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ABSTRACT 

The Minimum Mass Solar Nebula (MMSN) is a protoplanetary disk that contains the mini- 
mum amount of solids necessary to build the planets of the Solar System. Assuming that the giant 
planet s form ed in the compact configuration they have at the beginning of the "Nice model", 
Desch (2007) built a new MMSN. He finds a decretion disk, about ten times denser than the 
well-known Hayashi MMSN. The disk profile is almost stationary for about ten million years. 
However, a planet in a protoplanetary disk migrates. In a massive, long-lived disk, this question 
has to be addressed. With numerical simulations, we show that the four giant planets of the Solar 
System could not survive in this disk. In particular, Jupiter enters the type III, runaway regime, 
and falls into the Sun like a stone. Known planet-planet interaction mechanisms to prevent mi- 
gration, fail in this nebula, in contrast to the Hayashi MMSN. Planetary migration constrains 
the construction of a MMSN. We show how this should be done self-consistently. 

Subject headings: solar system: formation — accretion, accretion disks — methods: numerical 



1. Introduction 

The Minimum Mass Solar Nebula (hereafter 
MMSN) is the protoplanetary disk of solar com- 
position that has the amount of metals necessary 
to build the eight planets of the Solar System (and 
the asteroid belts). From the masses and composi- 
tions of the planets, a density of solids is derived at 
several locations of the disk. Then, the solar com- 
position is restored by adding gas, and a smooth 
protoplanetary disk density profile is derived. For 
a review of past works on th e MM SN concept, 
see the introduction of iDeschl ( 2007 ). The most 
famous version of the MM SN was provid ed by 
Weidenschillingl ()l977h and lHavashil (|l98ll ). As- 
suming that the giant planets have a rocky-icy core 
of about 15 Earth masses, and that Jupiter ac- 
creted all the solid material available between 1.55 
and 7 AU, Hayashi finds a surface density profile of 
E(r) = 1700( r/l AU)" 3 / 2 gx ri i" 2 , s imilar to the 
one given bv IWeidensch nhna (fl977h . Of course, 
the density profile obtained depends crucially on 
the position of the planets. Hayashi assumed that 
the planets formed where they presently orbit. 



A recent model explains several features of the 
Solar System (the Late Heavy Bombardment, the 
orbital distribution of the Trojans of Jupiter, the 
orbital elements of the giant planets. . . ) thanks to 
a late global instability in the o u ter Solar System 



dynamics dGomes et al. 2005; iMorbidelli et al 



2005: ITsiganis et al.ll2005l ). In this so-called "Nice 



model" , the four giant planets were in a compact 
configuration just after the solar nebula dissipa- 
tion. There fore, if one assumes th at th is model 
is true, the IWeidenschillina (|l977l ) and lHavashi 



(|198ll ) nebula i s out of date. In a recent, very nice 



article, Deschj ( 2007 ) constructed a new MMSN, 



assuming that the planets were formed in the disk 
at their starting position in the Nice model. He 
finds that the density distribution required to form 
the giant planets in this configuration, including 
the presence of an exterior disk of planetesimals, 
can be very well fitted by a steep power law density 
profile. In addition, he studies the time evolution 
of this disk. He finds a solution of a decretion 
disk, that can survive with an almost unchanged 
density profile in the giant planets region for a 
few million years. This enables the solid cores 
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of the four giants to reach their isolation masses, 
and then to slowly accrete their gaseous enve- 
lope. Consequently, this nebula is self consistent 
from the planetary formation point of view. The 
nice agreement between the initial position of the 
planets in the Nice model and a long-lived, almost 
power-law protoplanetary disk can be seen as a 
new plus point for the Nice model. 

However, planets in gas eous disks are subject to 
plane tary migration (see Papaloizou fe Terqueml 
120061 for a review). How to prevent Jupiter (and 
Saturn) from becoming a hot Jupiter (or Saturn) 
is a longstanding problem in Solar System for- 
mation, but some so l utions have been proposed. 
Masset fc Snellgrovd (|200lh have shown that a 



mean motion resonance between Jupiter and Sat- 
urn can enable them to resist the driving inwards 
by the disk, and to migr ate together outwards. 



Morbidelli fc Cridal (|200Tt> found disk parameters 



for which Jupiter and Saturn enter the 3:2 mean 
motion resonance and then hav e a negligible mi- 
gration. iMorbidelli et alj (|2007r ) then showed that 
Uranus and Neptune can also be saved in that 
case, by capture in resonance with Saturn. The 
four giant planets are in the end of the gas disk 
phase in a compact, fully resonant configuration. 
Two of these possible final configurations can lead 
to a late global instability of the dynamics, as re- 
quired in the Nice model. These works were based 
on a disk in which the density slope was low and 
the gas density at the loc ation of Jupite r was of 
the order of the one in the lHavashil (| 1 9 8 lh nebula. 

In this paper, we addre ss the q uestio n of the 
planetary migration in the iDeschl (|2007l ) nebula. 
In section [2j we describe our code, that enables 
a simulation of the entire disk, and of its global 
evolution. We also review the disk properties, and 
show how they should affect the planetary migra- 
tion. The results are presented in section [3] for 
a locally isothermal disk, and in section |4] for a 
disk where the energy equation is taken into ac- 
count. In every case, the planets are lost into 
the Sun. In section 15.1 1 w e try to rep roduce 
the IMorbidelli fc Cridal (|2007l ) result in the lDeschl 
torn disk, or in a colder disk with same density 
profile. In section [5.21 we focus on original results 
on Uranus and Neptune in some s imulations. For 
comparison, the migration in the lHavashil (|l98ll ) 
MMSN in presented in section [6l Some numeri- 
cal issues are discussed in section [7] Finally, our 



conclusions are presented in section and some 
prospectives are suggested. We show how plan- 
etary migration should be considered in the con- 
struction of a Minimum Mass Solar Nebula, and 
we draw some ideas towards a MMSN compati- 
ble with both planetary migration and the Nice 
model. 

2. Numerical settings, disk properties and 
expected migration 



2.1. Code description 

We use the cod e FARGO (iMassed teOOOallbh. in 
its 2D1D version^ (jCrida et al.ll2007t ). In this ver- 
sion, the classical 2D polar grid of the hydro-code 
(in which the planet-disk interactions are com- 
puted) is surrounded by a ID grid. In the ID 
grid, the disk is considered as axisymmetric, and 
only the radial component of the equations is com- 
puted ; this is computationally very cheap, there- 
fore the ID grid extension can be very large. The 
two grids are smoothly connected, so that the disk 
evolution is computed accurately over all the ID 
grid extension : the evolution of the inner and 
outer parts of the disk computed in the ID grid, 
influence the disk evolution in the 2D grid ; and 
reciprocally the planet-disk interactions computed 
in the 2D grid perturb the global disk evolution. 

This property of the code makes it the perfect 
tool to study the behavior of the giant planets in 
the MMSN. Indeed, giant planets generally open 
gaps in the disk, and then are locked in the global 
viscous evolution of the disk (type II migration). 
Most often, this drives the planets inwards. But in 
the Desch nebula, the disk is viscously spreading. 
Consequently, its global evolution should drive 
planets in type II mig ration outwards. However , 
Crida et al.1 (|2007l) and lCrida fc Morbidellil (|2007l) 
found that things are not that simple. If the gap 
is not completely empty, a corotation torque is ex- 
erted on the planet, gas can pass from the inner to 
the outer disk (or reciprocally) , and the planet can 
decouple from the disk evolution. Type II migra- 
tion can then proceed inwards in a decretion disk, 
or outwards in a viscous accreting disk. This can 
be seen only with the use of the 2D1D algorithm, 
that permits both an accurate 2D computation of 
the planet-disk interaction and the ID computa- 



x Now publicly available on http : //f argo . in2p 3 . f r] 
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tion of the global disk evolution. The global disk 
evolution also matters for the type III regime. 

The outer edge of the ID grid is loca ted at 
61 A U from the star, as prescribed by iDesch 
(|2007l ). The inner edge is located arbitrarily at 



0.2 AU. The inner and outer edges of the ID grid 
are open, allowing gas to flow out of the grid. 

The rings of the grid are logarithmically dis- 
tributed : 5r/r is constant through the two grids. 
The cells of the 2D grid are squared : Sr/r = 66. 

The 2D grid extends from 1 or 1.6 to 23.5 AU, 
covering the planets region. 

In the computation of the force of the disk on 
the planets, a part of their Hill sphere is excluded. 
To perform this, we calculate in the code the vec- 
tor force F p using the following expression : 



(1) 



cells 



where G is the gravitational constant, M p the 
mass of the planet, dS is the surface of the consid- 
ered cell, s is the vector from the planet to the cen- 
ter of the cell, s is its length, and s' = Vs 2 + e 2 is 
the smoothed distance to the planet. The smooth- 
ing length e is 0.6 th, where ry is the Hill radius 
of the planet. The term f(s) is our filter used to 
exclude the neighborhood of the planet given by : 



/(*) = 



exp -10 



0.6r H 



- 1 
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(2) 



It is a smooth increasing function from at s = 
to 1 when s — > oo, through 0.1 when s = 0.47 Th, 
1/2 when s = 0.6 rg, and 0.9 wh en s = 0.73 
(see Figure 2 of lCrida et al~ll2008l ). 



2.2. Disk and planets parameters 

The Solar Nebula found by iDeschl (|2007t ) has 



the following properties. Assuming that Jupiter 
formed at 5.45 AU, Saturn at 8.18 AU, Uranus at 
14.2 AU and Neptune at 11.5 AU, the gas density 
should be : 

E„(r) = 343(r/10AUT 2 - 168 g.cuT 2 (3) 
= 5.68 x 10~ 3 (r/lAUr 2 - 168 M©AU~ 2 . 
This is about 10 times more dense than the 



Havashil (| 19811) nebula at 5 AU, and 6.4 times more 
dense at 10 AU, and 3 times more dense at 30 AU. 



The steep density slope makes the disk a decretion 
disk, viscously spreading outwards, fed by the in- 
ternal parts. The steady state profile found by 
IDeschl (l2007h has a slightly different shape, but 
we use the power law for convenience. The two 
profiles are almost identical in the giant planets 
region. The outer edge of the disk is found to be 
located at 61 AU from the star. 

The temperature is 150(r/lAU)~ 429 K, which 
corresponds to an aspect ratio of 



(H/r) = 0.05(r/lAU) c 



(4) 



The viscosity is given by an a prescription 
(|Shakura fc Sunvaevl[l973h . with a = 4x 10~ 4 . 

With these characteristics, the disk profile re- 
mains almost unchanged for nearly 10 millions 
years, which leaves time for the solid cores of the 
planets to reach their isolation masses and accrete 
their atmosphere. However, this also leaves the 
planets time to migrate. 

In our simulations, the planets are initially lo- 
cat ed on c i rcular orbits at the position from which 
the IDeschl (|2007h disk is derived (see above). The 
planets masses are grown smoothly from to their 
present masses over 30 years at the beginning of 
the simulation. The planets are not accreting gas 
from the disk. 

During the first 100 orbits of Jupiter (1274 
years), the planets do not feel the disk potential, 
and therefore don't migrate. During this time, the 
planets launch a wake, open a gap, perturb the 
disk. This is for the disk to adapt to the planets 
potential and reach an equilibrium state. Then, 
the planets are released under the influence of the 
disk, and start their migration. 

2.3. Expected migration 

The aspect ratio is 8.1% at the location of 
Jupiter, and 9.1% at the location of Saturn. The 
Reynolds number 1Z — r 2 Q/v is 3.8 x 10 5 at the 
location of Jupiter, and 3.02 x 10 5 at the loca- 
tion of Saturn. Denoting V = | ^ + ^ , we 
find V = 0.876 + 0.131 w 1 for Jupiter, and 
V = 1.498 + 0.583 « 2 for Saturn. The unified 
criterion to open a gap of de pth 90% of the un - 
perturbed density is V < 1 (|Crida et al. 2006). 



Therefore, we expect Jupiter to open a non-empty 
gap, and Saturn not to perturb the density profile 
significantly. 
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Co nsequently the mechanism used bv lMorbidelli 
(120071 ) ;o prevent Jupiter and Saturn migration 
should not work : it requires that the two planets 
lie in a wide common gap in mean motion reso- 
nance. On the contrary, Saturn should be here in 
type I migration, like Uranus and Neptune. Given 
that these three planets are massive (15 to almost 
100 Earth masses), and that the disk density is 
high, their type I migration should be fast. 

The isothermal horseshoe drag exerts on the 
planet a torque proportional to the logarithmic 
gradient of the vortensity: (dln(E/i? )) / (dlnr) , 
where B is the second Oort constant |Wardlll992l 
Massed 12001 ) . In Kepler ian rotation, B = 0/4. 



fc Cvid& r/r = 0.01, so that it is divided in 270 rings, 



Thus, in the 



Deschl (|2007[ ) nebula, E/B oc r~ - 668 , 
and dln(E/B)/dlnr = -0.668. Due to the steep 
negative slope of the density profile, the corotation 
torque is negative, while it is zero if X oc r~ 1,5 
and positive if the slope is shallower. Such a 
disk density slope should therefore enhance the 
type I inwards migration of the planets with re- 
spect to more classical disks, in particular the 
Havashil (|l98lh MMSN. 

The high mass of the disk and the fact that 
Jupiter opens a small gap, makes it likely that 
Jupiter falls in the type III, runaw ay migration 
regime (jMasset fc Papaloizoul 120031 ) . Indeed, the 
mass of gas present in the coorbital region of 
Jupiter in the unperturbed disk is about 4 times 
the mass of Jupiter ; therefore, the coorbital mass 
deficit may exceed the mass of the planet, which 
is the condition for the runaway migration. 

In summary, in the Nebula found by iDeschl 
(|2007l ). the four giant planets of the Solar System 
should migrate inwards on a short timescale. To 
check this, and to study the possible interactions 
between the four planets in the nebula, we have 
performed numerical simulations, presented in the 
next sections. 

3. Models with locally isothermal equa- 
tion of state 

In this subsection, we use a locally isothermal 
equation of state P — c s 2 E, with c s the sound 
speed. The sound speed is a function of the dis- 
tance to the star given by c s 2 = (H/r) 2 GMQ/r cx 
T, which is not evolving with time. The energy 
equation is not computed. The 2D grid ranges 
between 1.6 and 23.5 AU, with a resolution of 



and 628 sectors. 

The planets are set according to section 12.21 
The migration path obtained is displayed in Fig- 
ure!]] As soon as released, Jupiter enters a type III 
migration regime and reaches 2 AU in 300 years. 
There, it stops due to interaction with the bound- 
ary of the 2D grid. If the grid had extended further 
inwards, there is no reason why Jupiter would not 
have reached the Sun. 

Saturn also migrates inwards, at first in type I 
migration because it hardly perturbs the density 
profile. It enters the type III migration regime 
only at t = 3600 years. Indeed, at this time, Sat- 
urn has reached 6.2 AU, where Psaturn = 1-88. 
Then, the dip dug by Saturn in the disk profile is 
deeper than before, the coorbital mass deficit can 
reach the mass of Saturn, and the runaway process 
can start. 

When Saturn reaches 3.7 AU, at t = 3900 
years, it is caught in the 5:3 mean motion reso- 
nance with Jupiter : for 300 years, the ratio be- 
tween their semi major axes is as/aj = (5/3) 2 / 3 
(where a denote the semi major axis, and the sub- 
script refers to the name of the planet). This 
breaks the runaway. At t = 4260 years, Saturn 
has a close encounter with Jupiter and is kicked 
out. Then, it starts again an inwards migration 
until it is blocked again by Jupiter : from 5600 
on, the ratio between their semi major axes is al- 
most constant as/aj 1.43. This may be due to 
a mean motion resonance or to indirect interac- 
tions, Jupiter perturbing the disk and thus the 
motion of Saturn ; in particular, the corotation 
torque should be positive and strong on the outer 
edge of the gap of Jupi ter, which repels Saturn 



[see iMasset et al.l l2006al) . The two planets then 
start to mi grate slightly outwards in a common 
gap, like in lMasset fc Snellgrovd (|200ll ). 

Neptune and Uranus are in the type I migra- 
tion regime. Between 2600 and 3600 years, Nep- 
tune is slowed down by Saturn : the ratio of the 
semi-major axis of the two planets remains con- 
stant equal to on/cls = 1-17 « (5/4) 2 / 3 during 
this period. Between 3060 and 3520, Uranus is 
also slowed down, with ajj /as — 1-39 « (5/3) 2 / 3 , 
and au/aN ~ (4/3) 2 / 3 . This shows likely cap- 
tures of Neptune in the 5:4 resonance with Sat- 
urn, and Uranus in the 5:3 with Saturn or the 4:3 
with Neptune. However, a mean motion resonance 
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can not be maintained if the involved bodies move 
away from each other. When Saturn accelerates 
its inwards migration, Neptune and Uranus stay 
behind. 

Once Saturn has moved inward, the two ice 
giants accelerate again their inwards migration. 
Neptune goes faster and reaches an other reso- 
nance with Saturn. From 4800 years on, their 
semi-major axis ratio is constant equal to (7/5) 2 ' 3 . 

At time 5900 years, Uranus reaches Neptune. 
The two planets then share the same average semi- 
major axis, equal to 1.5 2 / 3 as, which is character- 
istic of a 3:2 mean motion resonance with Saturn. 
Their configuration will be analyzed in more detail 
in section I5T21 

3.1. High resolution 



Because some authors (jD'Angelo et al.l 120051) 
have found that type III migration is resolution 
dependent, we have performed the same experi- 
ments with a better resolution of Sr/r = 10~ 25 . 
The 2D-grid is also extended inwards and ranges 
now from r = 1 AU to r = 23.5 AU. The 2D grid 
is now divided in 1000 rings and 1987 sectors. The 
ID grid is divided in 1811 rings between 0.2 and 
61 AU. 

The result is displayed in Figure [2] It is very 
similar to the previous case. Jupiter still migrates 
in type III migration as soon as released, and Sat- 
urn also undergoes a runaway migration episode, 
but a few hundreds of years earlier than before. 
The planets go further inwards due to the ex- 
tended 2D grid ; this shows that their stop is a 
numerical artifact, as expected. The global re- 
sult is not affected by the resolution. In fact, Fig- 
ures Q] and [2] look similar, and the same explana- 
tions hold. 

From 4100 years on, Saturn and Jupiter mi- 
grate slightly outwards in the 5:3 mean motion 
resonance, with as/aj — (5/3) 2 / 3 . 

Neptune and Uranus have the same migration 
rate in both resolution cases, before 3000 years. 
After 3000 years, they migrate faster in the high 
resolution case because Saturn is not present any- 
more to slow them down. At t = 5400, a close 
encounter happens between Uranus and Neptune 
(the distance between the two planets is smaller 
than 0.4 AU). Then, the two planets share the 
same semi major axis again. More detail on their 



configuration will be provided in section 15.21 

After 10 4 years, the four planets have escaped 
a fall into the Sun, thanks to the inner edge of the 
2D grid. They have reached a fully resonant con- 
figuration in which their inwards migration is pre- 
vented . This is roughly similar to lMorbidelli et aT 
(|2007l ). However, here the four giant planets lie 
within 7 AU from the Sun, with Jupiter being 
at 3.6 AU. This is incompatible with the present 
structure of the main asteroid belt. In addition, 
the last planet is not stopped by a resonance with 
Saturn or with the last but one planet, so that 
the two ice giants share the same orbit. This 
is because the disk is more massive here than in 



Morbidelli et al.l (|2007l ). and the gaseous torques 
are stronger than the more distant resonances. It 
is very unlikely that this final configuration may 
lead to a late global instability that would drive 
the planets on their present orbits. 

3.2. Getting rid of the hot Jupiter 

In the previous simulations, the migration of 
Jupiter was stopped because the giant planet 
reached the inner edge of the 2D grid. Conse- 
quently, what happens after this event is admit- 
tedly interesting, but somehow artificial. In a new 
simulation, the planets are removed when they 
reach a < 2 AU. 

The results are displayed in Figure [3] for the 
higher of two resolutions used above (Sr/r = 
10 -2 ' 5 ). Jupiter being not present anymore, Sat- 
urn goes on migrating inwards and finally reaches 
as < 2 AU as well. Then, Uranus and Neptune 
migrate inwards freely The differential Lindblad 
torque, responsible for type I migration, is pro- 
portional to Y^r^Vi 1 (H / r)~ 2 ; thus it increases as 
r -1 ' 74 when the planets approach the Sun. This 
explains their accelerating migration. Within 8000 
years, they have migrated inwards all the way to 
the grid edge too. The same happens at lower 
resolution. 

With this more realistic way of dealing with the 
planets at the inner edge of the 2D grid, the four 
giant planets are lost in less than 10 4 years. 

4. Non isothermal disk 

It has been shown recently that computing 
the energy equation can change dramatically the 
corotation torque, and thus the migration rate of 
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low mass planets (IPaardekooper fe Mellema 2006t 



Baruteau fc Masse? 2008 ). Planets in type II mi- 
gration should not be perturbed much, but this 
process may be criti cal for Neptune mass planets 
(<Klev fe Cridall2008l) . '. n the case considered in 
this paper, we expect this effect to be also signif- 
icant for Saturn and Jupiter because their gap is 
not empty. To check this, we have implemented 
in FARGO-2D1D the computation of the Energy 
equation as done in the FARGO- ADSG version. 

The Energy equation is the following, with e 
the surface density of internal energy of the gas, v 
the velocity vector, and P the pressure : 

De 

— + V(ev) = -PX7v + Q+ + Q- , (5) 
at 

where Q + and Q_ are heating and cooling terms, 
respectively. 

The heating term comes from the viscous heat- 
ing. The cooling term is given by a vertical black 
body emission : Q_ = 2er/?T 4 / kE, where gr is the 
Stefan-Boltzmann constant, T is the temperature, 
and k is the opacity. The opacity i s chosen such 
that the temperature profile given bv lDeschl (|2007t ) 
is an equilibrium profile of the unperturbed disk 
(in which Q + = 9E ^ 2 / 4 ) : 

8 

(6) 



8a R (H/r) c 
9 



with (H/r)o given by Eq. (gj). 

The cooling time r cool = e/Q_ is then (10/9a)fi~ 
T cooi = 2.78 x 10 3 Q _1 . For Uranus and Neptune, 
the horseshoe libration time Tub = 8irr p /(3£lx s ) 
i s 326 ft' 1 at 12 AU with x s = 1.16r py /q/(H/r) 
(|Masset et alj|2006bl Eq. (3)). The cooling time 



being large with respect to the libration time, the 
thermal effect on the corotation torque should sat- 
urate quickly and the migration of the ice giants 
should not proceed outwards like in lKlev fe Crida 
(l2008h . 



Our equation of state is, with 7 the adiabatic 
index (7 = 1.4) : 



p 



( 7 -l> 



(7) 



Now, the sound speed c s is not fixed but deter- 
mined at every time-step by the local temperature, 
which is evolving with time. 

The same experiments as in previous section are 
computed, with this energy equation and equa- 
tion of state. The result at lower resolution is 



shown in Figure |4j Jupiter and Saturn still un- 
dergo a type III migration episode. The migra- 
tion of Uranus and Neptune is slower than with 
a locally isothermal disk, but is still directed in- 
wards. In the end, the situation is basically un- 
changed with respect to previous section : Jupiter 
is blocked by the edge of the 2D grid, Saturn is 
blocked by Jupiter, Uranus and Neptune share a 
common orbit in resonance with Saturn. 

With a higher resolution (10 -25 ), and remov- 
ing the planets when they reach the inner edge 
of the 2D grid, one gets the migration paths dis- 
played in Figure [3 Jupiter and Saturn disappear 
below 1.3 AU within 4500 years. The migration 
of Uranus and Neptune is a bit slower than in the 
previous case. However, as soon as Saturn starts 
its type III inwards migration, the migration of 
the ice giants accelerates. Then the migration of 
Neptune is similar to the previous case (with a 
800 years delay), while that of Uranus is much 
slower (au w 2 X 10" 4 AU.year" 1 ). After 8800 
years, Neptune reaches 1.1 AU and is declared 
lost. Uranus goes on migrating slowly inwards, ac- 
celerates after Neptune is suppressed, and reaches 
1.3 AU after 15200 years. 

5. Interactions between the planets 
5.1. Jupiter and Saturn 

Previous sections have shown that if the 4 gi- 
ant pl a net o f the Solar System are placed in the 
Deschl (|2007l) nebula in the starting configuration 
of the Nice model and released simultaneously, 
they migrate inwards and disappear in the Sun 
in about 10 4 years. In this section, we try 
to prevent this by releasing the planets in se- 
quence, and varying the disk aspect ratio. We 
aim at reproducing the best candidate mecha- 
nism so far to preve nt the inwards migra t ion o f 
Jupi ter and Saturn : Masset fe Snellgrovel ( 2001 ) 
and iMorbidelli fe Cridal (|2007fT have shown that 
when Jupiter and Saturn orbit in a mean motion 
resonance in a common gap, they can avoid type II 
migration or even migrate outwards. 

Therefore, we have performed simulations in 
which Jupiter is held on a fixed circular orbit 
longer than Saturn and the ice giants. The three 
lighter planets migrate inwards. Once Saturn's 
semi major axis has reached a constant value in 
terms of Jupiter semi major axis (which betrays 
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a resonance), we release Jupiter as well. Immedi- 
ately, Jupiter migrates inwards in type III migra- 
tion and disappears ; and the three other planets 
migrate inwards again. This is because the con- 
sidered disk is too thick to enable Jupiter to open 
a deep and wide gap, embracing also Saturn. 

So, we decrease the aspect ratio to H/r = 0.04. 
As the viscosity is given by an a-prescription, the 
viscosity is also smaller. In fact, "Pjupitor = 0.467. 
The planets are maintained on a fixed orbit for 
5250 years (500 Jupiter orbits). At this time, 
Jupiter has opened a deep and wide gap, the den- 
sity at the bottom of which is a bit less than 2% 
of the unperturbed density. Saturn has opened 
a partial gap, of depth a half of the unperturbed 
density, and the two gaps have merged (see the 
solid profile in Figure [51 compared to the dot- 
dashed initial profile ) . This kind of gap suits the 
Masset fe Snellgrovel ( 200lh mechanism: Jupiter 
and Saturn lie together between the inner and the 
outer disk. 

First, only Saturn is released. One expects than 
Saturn, repelled by the outer disk, would migrate 
inwards until it encounters a mean motion reso- 
nance with Jupiter. However, Saturn goes imme- 
diately in type III migration outwards and reaches 
10 AU. Twenty orbits later, Jupiter, Uranus and 
Neptune are released as well. Jupiter runs away in- 
wards and reaches 2 AU in about a hundred years. 
If the 4 planets are released simultaneously after 
500 orbits of Jupiter, the same thing happens. 

In an other attempt, Uranus and Neptune are 
not considered, and the aspect ratio is set con- 
stant equal to 0. 05, like in the almost stat ionary 
solution found bv lMorbidelli fe Cridal (|2007h . Sat- 
urn is placed on a circular orbit at 10 AU, and let 
free to migrate after 2550 years (200 Jupiter or- 
bits), while Jupiter is still held on a fixed orbit 
at 5.45 AU. Saturn migrates inwards in type III 
migration, and then is blocked by Jupiter. Then, 
we release Jupiter as well after 5100 years (400 
Jupiter orbits). Again, Jupiter runs in type III 
migration into the Sun as soon as released. 

Decreasing even more the viscosity and aspect 
ratio of the disk may help, but H/r — 0.04 and 
a = 4 x 10~ 4 are already rather small com- 
pared to the standard values inferred from obser- 
vations. In fact, the r e are two reasons why the 
Masset fc Snellgrovel ( 2001 ) mechanism can not 



that it is cool enough to enable Jupiter to open a 
wide gap : 

1. The density of this disk is so high that 
Jupiter (or Saturn) can easily enter the run- 
away type III migration regime. In that case, 
Jupiter is driven inwards faster than Saturn, 
and the putative resonance is broken. 



work in the considered disk, even if one assumes 



This d isk is a decretion disk. The lMasset fc Snellgrove 
([20011 ) mechanism prevents a pair of planets 
from inwards type II migration. It works 
when the innermost planet is more mas- 
sive than the outer one, and repels the later 
outwards against the accreting outer disk. 
Here, the global evolution of the disk is di- 
rected outwards. A more massive Saturn 
than Jupiter would be needed to prevent 
the outwards type II migration of the two 
planets. 



5.2. Uranus and Neptune 

In the above simulations, when the migration 
of Jupiter is artificially stopped at the edge of the 
2D grid, Uranus and Neptune always end on the 
same orbit, captured in the same mean motion 
resonance with Saturn. This strange configuration 
is worth a more detailed analysis. 

In the first case presented (locally isothermal 
disk, with resolution Sr/r — 10 -2 ), the distance 
between the two planets - shown as crosses in Fig- 
ure 0- is bounded by 0.1 AU after t = 6000 
years. The Hill radius of the planets is r# = 
a p (M p /3) 1 / 3 = 0.14 AU at a p = 5.5 AU. This in- 
dicates a satellite motion of the planets. This kind 
of configuration requires a lot of damping because 
the planets were of course unbound gravitationally 
at the beginning. The damping is here provided 
by the disk, and is visible between 5900 and 6000 
years in Figure[7]: the maximum distance between 
the planets decreases until it becomes smaller than 
the Hill radius. The high density of this nebula 
makes such a satellite capture possible. 

In the second case (locally isothermal, 5r/r = 
10~ 2 ' 5 ), the planets are not satellites of each other. 
Here, the distance between Uranus and Neptune is 
slightly oscillating around 2 AU after 6000 years. 
They share the same orbit, with an eccentric- 
ity smaller than 0.04 and Uranus leading about 
30° before Neptune : the angle between Uranus 
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and Neptune, seen from the Sun, is plotted as 
crosses in Figure El with respect to the right y- 
axis. The two planets are not at a Lagrange point 
of the other one. This unusual configuration may 
be due to the mean motion resonance with Sat- 
urn, or more probably to the strong dissipation by 
the disk, that modifies the local dynamics. Once 
the disk will have evaporated, these two planets 
should have a tadpole motion around one of their 
L4/L5 Lagrange points. Such a configur ation have 
alread y been studied , for in s tance by lThomin.es 
(|2005[) ; iBeauee et all (|2007t) ; ICresswell fc Nelson 
(|2008l ). Generally, it is assumed that a giant 
planet captures a migrating terrestrial planet, that 
may then grow. In contrast, here the effect of the 
giant planet is to put two 15 Earth mass planets 
together on the same orbit. This is a new way of 
forming Trojan planets. 

In the third case (with energy equation, and 
Sr/r — 10 -2 ), the distance between the two ice 
giants (shown as the + symbols in Figure [4j set- 
tles to about 2 times their semi-major axis. This 
indicates that the two planets are in opposition, 
librating around their L3 point. Once again, this 
configuration is not possible in the absence of gas 
dissipation or of the other planets, because the L3 
point is unstable. If the other planets and the disk 
disappear, these two planets should have coorbital 
horseshoe orbits. 

These three original outcomes open new pos- 
sibilities for extra-solar systems. They do not fit 
the Solar System configuration, though. Whether 
or not a scenario like the Nice model may apply 
from such 1:1 resonances, is not clear. In particu- 
lar, the configurations with Uranus and Neptune 
at opposition or at 30 degrees on the same orbit 
are probably unstable on the short term in absence 
of the gas disk. This is thus unlikely to yield to 
a late instability, as required. In the case where 
the two planets are satellite of each other at a dis- 
tance d = 0.1 AU, the energy required to separate 
them, (GM v M N /d - GMuM N /r H )/2, is roughly 
equivalent to the energy required to eject a Moon 
size body out of the Solar System from a 10 AU 
orbit. Therefore, it is likely that they stay in this 
satellite configuration. 



Havashil ()198ll ) nebula 



6. Migration in thelWeidensch illing (1977^) 



For comparison, the same experiment as in sec- 
tion [3] has been performed in the standard Min- 
imum Mass Solar Nebula. The initial density is 
the one quoted in introduction, that is : 

S(r) = 1.9125 x 10~ 4 (r/l AUT 3/2 M .AU~ 2 . 

(8) 

This profile appears in Figure [6] as the straight, 
gray, dotted line. 

The equation of state is locally isothermal. The 
aspect ratio is given by Eq. (j4]), and the viscosity 
is given by an a-prescription with a — 4 x 10 -4 . 
The planets are initially placed on circular or- 
bits at their present position: aj — 5.2 AU, 
a s = 9.6 AU, au = 19.2 AU, = 30 AU. The 
resolution is 10 -2 , the 2D grid extends radially 
from 2.15 to 52 AU, while the ID grid covers the 
range [0.312 ; 104] AU. 

The planets do not feel the disk potential dur- 
ing 400 Jupiter orbits (4740 years), and are then 
released. The evolution of the semi major axes of 
Jupiter and Saturn are displayed in Figure [5) In 
spite of the fact that the shape of the gaps opened 
by the planets are similar to the ones in the lDesch 
(2007) disc (the wi dth and depth dep end only on 



M p , H/r and v, see lCrida et al.ll2006l ). no type III 



migration take place, because this disc is not mas- 
sive enough. Saturn migrates faster than Jupiter ; 
about 2600 years after the release, the two planets 
enter in 2:1 mean motion resonance, which slows 
Saturn down. After 25000 years, the migration 
stops at aj = 4.75, a s = 7.57 AU. 

In an other, longer term simulation, we change 
the aspect ratio to H/r = 0.05 and the viscos- 
ity to a = 4 x 10 -3 . The low aspect ratio en- 
ables Jupiter to open a significant gap, and Sat- 
urn a shallow one. The density profile after 400 
Jupiter orbits can be seen as the bottom, gray, 
dashed line in Figure [6) The migration paths are 
displayed in Figure [9] for 10 5 years. Jupiter mi- 
grates inwards in type II migration, while Sat- 
urn migrates faster. After 6400 years, Saturn 
encounters the 2:1 mean motion resonance with 
Jupiter, but is not caught in it. At about 10 4 
years, Saturn encounters the 5:3 mean motion res- 
onance with Jupiter; from this time on, as/aj 
(5/3) 2 / 3 . The gaps of the two pla nets merge. 
Then, the Masset k, Snellgrovei ( 200ll ) mechanism 
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takes place, and the two planets have very li t- 
tle migration, like in iMorbidelli fc Cr ida| (l2007h : 
their semi major axes don't change by more than 
1 AU in 10 5 years. 

Note that the stop of the migration of Jupiter 
at 4.5 (or 4.75) AU is not due to the edge of the 2D 
grid : the radius of the inner edge of the 2D grid is 
about a half of the semi major axis of Jupiter. The 
presence of Saturn is responsible for the salvation 
of Jupiter. To check this, the same simulation 
has been performed with Jupiter alone, without 
the other three planets. The migration path is 
displayed in Figure [9] as the orange dashed line. 
In this case, Jupiter has an unperturbed type II 
migration, and reaches 3 AU in 40 000 years, 2.7 
AU in 50 000 years. 

For comparison, the migration path of Jupiter 
in Figure[5]is displayed also in Figure[9]as a dotted 
red line. It looks like a free fall towards the Sun. 

Uranus and Neptune are migrating inwards in 
type I migration. Their migration speed is slower 
than in the previous cases, because the density 
is moderate and the density slope is such that 
the isothermal corotation torque is zero. When 
approaching Saturn, Uranus slows down and al- 
most stops between 9 and 10 AU. Then, Neptune 
is caught in the 2:3 mean motion resonance with 
Uranus at 8.5 x 10 4 years : from this date on, their 
semi-major axis ratio is 1.5 2 / 3 . 

6.1. Conclusion 



Discussion 



7.1. Energy Equation and cooling law 



Migration also occurs in the Havashil (1981) 
nebula, and it is not possible that the giant 
planets stay where they form unt il the present 
time. Therefore, the lHavashil (1981) nebula is not 
more consiste nt with planetary migration than the 
Deschl (|2007l ) one. However, no runaway migra- 
tion takes place, and the Jupiter-Saturn pair can 
avoid migration all the way to the Sun for some 
reasonable disk param eters. Then, a sc e nario like 
the one presented in IMorbidelli et al. (|2007l) is 
possible. In this article, the authors find 6 fully 
resonant configurations of the four giant planets 
in which the migration is prevented. After disk 
dissipation, two of these configurations can lead to 
a late global instability of the planetary system, 
compatible with the Nice model. 



To slow down migration in the iDeschl ((20071) 
nebula, one can think of taking the energy equa- 
tion into account with an other opacity law. This 
would change the thermal structure of the disk. 
With a shorter cooling time, the temperature and 
the aspect ratio decrease. This has basically two 
effects : (i) Saturn and Jupiter would open a 
deeper gap, but our locally isothermal simulations 
with H/r ^ 0.05 show that it may not be enough 
to prevent their migration ; (ii) the thermal horse- 
shoe effect would not saturate if t coo i ~ rub , which 
could lead to a positive corotation torque on the 
planets, particularly Uranus and Neptune. 

This should be the subject of a future, dedi- 
cated study. However, it is rather unlikely that 
the four planets stay exactly on their initial or- 
bit until the disk dissipates, in order for the Nice 
model to take place. 

7.2. Accretion 

In all the simulations presented so far, the plan- 
ets were not allowed to accrete. This may appear 
unrealistic at first sight, but on the other hand, 
accretion must have stopped anyhow before the 
masses of the planets exceed the present masses. 
This is certainly a critical issue for the forma- 
tion of Jupiter and Saturn, because the accretion 
timescale for so massive planets should be very 
short. In both density profiles studied here, free 
accretion of the planets leads to too high masses. 
This open problem requires a complete, dedicated 
study, that is beyond the scope of this paper. We 
make the assumption that some mechanism pre- 
vented accretion to go on when Jupiter and Sat- 
urn reached their actual masses. However, a few 
comments are in order. 

In the simulation presented in Figure [H the gas 
density at the bottom of the gap created by Sat- 
urn is of the order of 3 times the one at the bot- 
tom of the gap of Jupiter. If this situation lasts a 
few hundred thousand years, Saturn may accrete 
more gas than Jupiter on t he long term. This 
would be a problem for the iMasset fc Snellgrove 
(|200lh mechanism: a lighter outer planet than 
the inner one is required for the migration to be 
slowed down, stopped, or reversed. In the Solar 
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System though, we know that the mass of Sat- 
urn somehow remained smaller than a third of the 
mass of Jupiter. This mass ratio allows the mi- 
gration of the pair to be negligible over the disc 
life time, with suitable disk parame ters. There- 
fore, the Ma sset fc Snellgrovd ( 2001 ) mechanism 
applies. 

An other issue with accretion c oncer ns type III 



migration. iD'Angelo k. Lubow (|2008h have ob 



served that allowing the planet to accrete, may 
prevent the runaway migration. However, this is 
most likely because their accreting planets grow 
too fast from low mass planets in type I migration 
to high mass planets in type II migration in a clean 
gap. With accretion turned on, but the planet 
mass artificially kept constant to 3 x 10~ 4 M„, 
they find that type III migration occurs, which 
shows that accretion should not perturb too much 
the corotation torques responsible for the run- 
away. To check this, we have restarted the sim- 
ulation presented in section 13.11 and Figure [2] at 
the time where the planets are released, turning 
accretion by Jupiter and Sat urn on. Acc retion is 
computed using the recipe bv lKlevi l|l99flh . with a 
timescale to accrete all the gas within 0.4577/ of 
the planet equal to 1.6 year, that is 1/8 Jupiter 
orbit. The evolution of the mass and semi ma- 
jor axis of Jupiter are displayed in Figure [10] as 
the curves with + symbols, compared to the solid 
line taken from Figure [2] Type III migration is 
confirmed (and even enhanced), even if Jupiter 
reaches 5 Jupiter masses in a century. 

For all these reasons, and even if a deeper study 
of the accretion processes would certainly be valu- 
able, we think that our results are robust as far as 
this issue is concerned. 

8. Summary and conclusion 



In this paper, we have confronted the iDesch 
(|2007l ) Solar Nebula to the planetary migration 
issue. This disk has a high gas density, a high 
density slope, and this profile remains almost un- 
changed for 10 million years. Therefore, the ques- 
tion of the planetary migration should not be 
eluded. The 2D1D version of FARGO is partic- 
ularly suited to such a study. The physical pa- 
rameters of the disk are fixed, but a broad range 
of numerical parameters have been studied (res- 
olution, energy equation, fate of the planets that 



reach the inner edge of the 2D grid), as well as the 
influence of the aspect ratio. 

The high density slope gives a negative non 
thermal horseshoe torque, that accelerates even 
more the inwards type I migration of Uranus and 
Neptune. The high density permits type III mi- 
gration of massive planets. In particular, we can 
not avoid an extremely fast inwards migration of 
Jupiter. In spite of several attempts, it seems to be 
impossible to prevent a loss of all the giant planets 
in less than 2 x 10 4 years, which is 2 to 3 orders of 
magnitude shorter than the disk life time. Taking 
into account the energy equation with a cooling 
compatible with the assumed temperature struc- 
ture of the disk does not change the result. 

Preventing the migration of the giant planets 
requires interactions between the planets. How- 
ever, as this disk is a mass ive, decretion disk, the 
Masset fc Snellgrovel ( 2001 ) mechanism to prevent 



Jupiter and Saturn from migrating inwards can 
not be applied. Shortly said, we conclude that 
the lDeschl (|2007t) nebula is incompatible with our 
present knowledge of planetary migration. 

So, we are facing a problem. Either this new 
MMSN makes the planetary migration in the So- 
lar System an even more critical issue and makes 
obsolete the former sol utions to prevent migra- 
tion of the giant planet s ([Morbidelli fc Cr ida 200l 



Morbidelli et alJ 120071 ). or the existence of plane- 



tary migration makes this MMSN questionable (in 
particular the apparently unavoidable type III mi- 
gration of Jupiter). 

A possible solution has been laid out in sec- 
tion [6l If one believes that the Nice model is 
true, it is clear that the outer planets should have 
ended the disk phase in a more compact config- 
uration than the present one. However, the mi- 
gration of the planets can account for a formation 
on a larger radial range, followed by a compaction 
of the configuration of the giant planets, like in 
Figure [HI Then, a fully resonant configuration can 
be achieved, th at can be compat i ble wi th the Nice 
model, like in iMorbidelli et alJ (|2007h . In addi- 
tion, the solid material that built the giant planets 
may come not only from the region around their 
respective orbit s : dust drifts inwards in a proto- 
planetary disk ( Weidenschillind 1977 ). and small 
bodies migrate as well, so that the giant planets 
region may be replenished in solids by the outer 
parts of the disk. This enables the formation of 



10 



Jupiter, Saturn, Uranus and Neptune in a less 
massive disk than the Desch nebula, to avoid the 
type III migration of Jupiter. In se ction[6l we have 
seen that this is possible in the old lHavashil (1981) 
nebula. However, there is no reason why the gi- 
ant planets should have formed exactly where they 
presently orbit ; in particular the presence of an 
outer cold disk of planetesimals (required in the 
Nice model) is problematic if Neptune formed be- 
yond 25 AU. Thus, a new construction of a MMSN 
is needed, that takes into account the Nice model 
and the plan etary formation constraints, like in 
Deschl (|2007n . and also the migration of planets 
and planetesimals. 

In any case, our results show that planetary mi- 
gration should be considered in the construction 
of a MMSN. The location where the planets form 
determines the gas density profile of the nebula, 
which determines the migration path of the plan- 
ets, which drives the planets to a new position 
after the disk dissipation. This final configura- 
tion, and not the initial one, should be compatible 
with the Nice model (or with the present configu- 
ration if one does not believe that the Nice model 
is true). This idea requires a detailed study, that 
is beyond the scope of this paper, but the results 
presented here advocate for such a self-consistent 
construction of the Minimum Mass Solar Nebula. 
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Fig. 1. — Migration path of Jup it er, S aturn, 
Uranus and Neptune in the iDeschl (|2007l ) neb- 
ula, with locally isothermal EOS and resolution 
of 10- 2 . 




Fig. 



Migration path of Jup it er, S aturn, 



Uranus and Neptune in the IDeschl (|2OO70 neb- 
ula, with locally isothermal EOS and resolution 
of 10~ 2,5 (curves, left y axis). The crosses show 
the angle between Uranus and Neptune when they 
share the same orbit (right y-axis). 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 
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Fig. 3. — Migration path of Jup it er, S aturn, 
Uranus and Neptune in the iDesch ( 2007 ) neb- 
ula, with locally isothermal EOS and resolution 
of 10~ 2 ' 5 . The planets are removed from the sim- 
ulation when they reach the 2D-grid inner edge. 



2000 4000 



6000 8000 1 0000 
time [years] 



12000 14000 16000 



Fig. 5. — Migration pat h of Jupiter, Saturn, 
Uranus and Neptune in the iDeschl (|2007l) nebula, 
with Energy equation computed, and resolution 
of 10~ 2 ' 5 . The planets are removed when they en- 
counter the inner edge of the 2D grid. 
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time [years] 

Fig. 4. — Migration pat h of Jupiter, Saturn, 
Uranus and Neptune in the IDeschl (|2007l ) nebula, 
with Energy equation computed, and resolution 
of 10~ 2 (curves). Crosses: Uranus- Neptune dis- 
tance. 
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Fig. 6. — Top black lines: Initial density profile 
(dot-dashed, straight line) and pr ofile af t er 500 
Jupiter orbits (solid curve) in the IDeschl ( 2007 ) 
nebula with H/r = 0.04 (see section l5~Tj) . Bot- 
tom gray lines : Initial density profile (dotted, 
straight line) and pro file after 400 J upiter orbits 
(dashed curve) in the lHavashil (|l98ll ) nebula (see 
section [5]) . 
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Fig. 7. — Migration path of Uranus and Nep- 
tune taken from Figure [1] (two top curves labeled 
"auranus" and "a No ptunc" respectively). The size 
of Neptune Hill radius is the bottom light blue 
curve ("ru (Neptune)" ). The distance between 
the two planets is shown with + symbols, labeled 
"dist U-N"). 
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Fig. 9. — Migration path of Jupiter, Saturn, 
Uranu s and Neptune in a disk with the lHavashi 
(|l98lh density profile, H/r = 0.05, a = 4 x 10" 3 . 



The Equation of State is locally isothermal. The 
resolution is Sr/r = 10~ 2 . The orange dashed line 
corresponds to a simulation where Jupiter alone is 
migrating in the same disk. The red dotted curve 
is taken from Figure El for comparison. 
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Fig. 8. — Migration path of Jupiter (bottom red 
curve) a nd Satur n (top green curve) in a disk 
with the lHayashil (1981) density profile, H/r = 
0.05(r / lAU) °- 2855 and a = 4 x 10~ 4 like in the 
iDeschl (|2007t ) nebula. The Equation of State is 
locally isothermal. The resolution is Sr/r — 10~ 2 . 
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Fig. 10. — Bo ttom panel : Migration path of 
Jupiter in the IDeschl ( 2007 ) disk, with locally 
isothermal equation of state and resolution 10 -2 5 
without accretion (solid line, taken from Figure[2]), 
or with accretion with timescale 1.6 year (solid line 
with + symbols). Top panel : evolution of mass of 
the planet in the accreting case. 
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